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Cells modify their volume in response to changes in osmotic pressure but it is usually assumed that other ac¬ 
tive shape variations do not involve significant volume fluctuations. Here we report experiments demonstrating 
that water transport in and out of the cell is needed for the formation of blebs, commonly observed protrusions 
in the plasma membrane driven by cortex contraction. We develop and simulate a model of fluid mediated 
membrane-cortex deformations and show that a permeable membrane is necessary for bleb formation which 
is otherwise impaired. Taken together our experimental and theoretical results emphasize the subtle balance 
between hydrodynamics and elasticity in actively driven cell morphological changes. 


Cells can change their shape to explore their environment, 
communicate with other cells and self-propel. These macro¬ 
scopic changes are driven by the coordinated action of lo¬ 
calized motors transforming chemical energy into motion. 
Active processes in biological systems can be linked to a 
large variety of collective non-equilibrium phenomena such 
as phase-transitions, unconventional fluctuations, oscillations 
and pattern formation oHa. A vivid example of actively 
driven non-equilibrium shape fluctuations is provided by cel¬ 
lular blebs, the rounded membrane protrusions formed by the 
separation of the plasma membrane from the cortex as a result 
of acto-myosin contraction (Ml. 

Blebs occur in various physiological conditions EUSl, as 
for instance during zebraflsh embryogenesis ITtUTTII. or cancer 
invasion 0. While some questions concerning the mecha¬ 
nisms governing bleb formation and its relation to migration 
have been resolved dClIEHni, key aspects of bleb mechan¬ 
ics remain unclear. Geometrical constraints dictate that ac¬ 
tive shape changes associated with blebs should necessarily 
involve either fluctuations in the membrane surface or in cel¬ 
lular volume, and possibly both. It is generally believed, how¬ 
ever, that the cellular volume is not signiflcantly altered during 
bleb formation, so that the cell is usually considered incom¬ 
pressible iiiiiiiii. Yet, experimental evidence in vitro sug¬ 
gests that aquaporins (AQPs), a family of transmembrane wa¬ 
ter channel proteins ca, are involved in cell migration CH- 
ira and blebbing (T^EOl. The implied significance of fluid 
transport through the membrane suggests that an interplay be¬ 
tween hydrodynamic flow and active mechanics has an impor¬ 
tant but still unclear role in blebbing. In this letter, we reveal 
the role of the membrane permeability in the formation, ex¬ 
pansion and retraction of cellular blebs. We show by direct 
experiments in vivo and numerical simulations that bleb for¬ 


mation involves volume fluctuations, considerable water flow 
through the membrane, and relatively smaller surface fluctua¬ 
tions. 

Experiment: One of the limitations impeding the experi¬ 
mental studies of the bleb dynamics is the lack of proper tools 
to generate high-resolution spatial-temporal data of bleb dy¬ 
namics. This is due to the fact that the time scale of bleb for¬ 
mation is relatively short (about 1 minute starting from initia¬ 
tion of the bleb to its retraction), which requires fast imaging 
and photostable markers. Here, we create an improved mem¬ 
brane marker ED which we inject in one-cell stage zebraflsh 
embryos (see Supplemental Material for experimental meth¬ 
ods 1221 ). Together with wild type (WT) zebraflsh PGCs, a 
well studied biological model to investigate blebbing in vivo 
ITtUTTI. we also consider cells expressing dominant-negative 
Rho Kinase mutant (DN-ROK) which inhibits acto-myosin 
contractility suppressing blebbing activity 0. Zebraflsh ex¬ 
presses a large number AQPs contributing to the water perme¬ 
ability of the membrane 1231 . Here we focus on AQPl and 
AQP3, the most ubiquitously expressed aquaporins 12211231 . 
To asses their role in volume change and blebbing, we con¬ 
sider PGCs with AQPl and AQP3 overexpression (AQP-f) and 
knockdown (AQP-). 

The imaging of blebbing is done during 12-16 hours post 
fertilization, when we record time-series of confocal images 
for a large number of cells. Sequences of image stacks are 
then processed using the 3D Active Meshes algorithm im¬ 
plemented in the Icy software 1241 . The algorithm performs 
three-dimensional segmentation and tracking using a trian¬ 
gular mesh that is optimized using the original signal as a 
target. Prom the resulting three dimensional mesh one can 
then measure the cell volume and its surface area (see also 
(23). In Fig. we illustrate representative phenotypes 
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of PGCs under different conditions (see also the correspond¬ 
ing movies (221). These observation show that WT cells dis¬ 
play a marked blebbing activity that, as expected, is strongly 
suppressed when active contraction is hindered, as in DN- 
ROK cells il. Remarkably, we also observe a strong re¬ 
duction in bleb activity in AQP- cells, where water flow is 
hindered. In contrast, water flow enhances blebbing as man¬ 
ifested by the presence of larger blebs in the AQP-f condi¬ 
tion. To quantify these qualitative observations, we measure 
the cell volume V and its surface S sampling the results over 
a large number of time-frames taken on different cells. To ac¬ 
count for cell-to-cell variability, we consider relative volumes 
AV/V = (V -_V)/V and surfaces AS/S = (S - S)/S 
changes, where V (H) is the time-averaged volume (surface) 
of each cell. The average value of the volume V does not 
change signiflcantly for the four cases (2 ^ . 


between AQP- and DN-ROK are not 1^ . Relative surface 
fluctuations are small that volume fluctuations in the WT case 
and are further reduced for DN-ROK and AQP- and slightly 
increased for AQP-f (see Fig. [T]:). We also checked that rela¬ 
tive surface and volume fluctuations are correlated (221, sug¬ 
gesting a direct link between blebbing activity and volume 
fluctuations induced by water transport. Suppressing water 
flow has the same effect as suppressing active contraction, in 
both cases blebs are hindered. Furthermore, the volume fluc¬ 
tuations we observe follow closely bleb expansion and retrac¬ 
tion as illustrated in Fig. Expansion of a bleb correspond 
to a visible volume increase while when a bleb retracts the 
volume decreases. Here we concentrate on AQP-i- PGC since 
the blebs are distinctly visible and the analysis clearer, but the 
same result holds for WT cells where, however, several blebs 
may form and retract simultaneously. 



FIG. 1: (Color online) Volume and surface fluctuations during bleb¬ 
bing are controlled by AQPs. a) Representative phenotypes for bleb 
formation in PGCs under different conditions: wild type (WT), cells 
expressing DN-ROK mutant which impairs contractility by interfer¬ 
ing with acto-myosin contraction (DN-ROK), cells where AQPla 
and 3a is suppressed (AQP-), cells with over-expression of AQPla/3a 
(AQP+). The scale bar corresponds to lO/xm. The cumulative dis¬ 
tribution of relative volume (b) and surface (c) fluctuations for the 
four conditions illustrated in a) together with a Gaussian fit (dashed 
lines). The distributions are sampled over N different time-frames 
corresponding to n different cells. (WT: n = 15 cells and N = 800 
time-frames; DN-ROK: n — 2 and N — 155, AQP-: n — A and 
N = 180; AQP+: n = 7 and A = 183). 

In Fig. [^, we report the cumulative distribution of rela¬ 
tive volume changes which indicates significant fluctuations, 
reaching up to 10%, in the WT case. Volume fluctuations 
are strongly reduced for the DN-ROK and AQP- cases, while 
they are enhanced in the AQP-i- case. The relative volume and 
surface distributions themselves are well described by Gaus¬ 
sian statistics, as also confirmed by a Kolmogorov-Smirnov 
test (22l. A statistical test also indicates that the differences 
between WT and both AQP- and DN-ROK are significant 
(p < 0.01), but the differences between WT and AQP-f and 


FIG. 2: (Color online) Bleb formation is directly related to volume 
changes, a) Time lapse of a single AQP+ cell. Two dimensional 
sections are shown in the upper panels and the corresponding recon¬ 
structed three-dimensional meshes in the lower ones. Arrows indi¬ 
cate blebs. The scale bar corresponds to lOpm. b) Evolution of the 
volume for the same cell. Time points corresponding to the panels 
in a) are denoted in yellow. An increase in volume is observed in 
correspondence with each newly formed bleb. 

Model: In order to better understand the physical role of 
water flow in bleb formation, we resort to numerical sim¬ 
ulations of a two-dimensional model of the biomechanics 
of cortex-membrane deformations including fluid transport 
through the plasma membrane. Several existing computa¬ 
tional models for cellular blebs simulate the detachment and 
expansion of the membrane due to the active contraction of 
the cortex, assuming cell volume conservation (25^281 . Here, 
we relax this assumption by introducing and varying the mem¬ 
brane permeability a (291 in a model based on the immersed 
boundary method (3Ql in the Stokes approximation (311 . con¬ 
sidering a contracting discretized elastic cortex coupled to an 
elastic membrane. 

We describe both the membrane and the cortex by a set 
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of discrete nodes connected by springs on one dimensional 
closed curves parametrized by their initial arc length s with 
s G [0,27r]. Using Lagrangian coordinates, the position of 
node i is denoted by . (t) (for the membrane) and Vc-{t) 
(for the cortex). Initial positions are chosen to be r^. = 
{vm cos(s), sin(s)) and = {vc cos(5), Tc sin(s)) where 
Vm and Vc are the membrane and cortex radii, respectively. 
The interactions between the nodes on their respective curves 
include nearest-neighbor (NN) and three-body interaction 
terms to account for stretching and bending energies. We 
also model the cortex-membrane interface by a set of springs 
with random stiffness k^c (Nm“^), drawn from a uniform 
distribution. Disorder in the stiffness represents at a coarse¬ 
grained scale the random arrangement of cortex-membrane 
linker proteins. The energy of the system can be decomposed 
SiS (o = + ^int • The first two terms are given by 


2e2 




) - cos(6>°^_^ . 


+ 
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( 1 ) 


where x corresponds to either membrane (m) or cortex (c), 

I is the distance between node i and i + 1, Cx is the equi¬ 
librium distance between NN nodes, (N/m) is the stiffness 
coefficient, (J) is the bending coefficient, and 

are the angles between the triplets (i — 1, i, i -h 1) 
in deformed and equilibrium configurations, respectively. The 
interaction energy is given by 4„t = X), ^mc(|rmci I - 0^/2. 
where |rmcj and I are the distances between the membrane 
and cortex node with the same index i in the deformed and 
equilibrium configurations. A non-vanishing rest length Cc for 
the cortex element is needed to prevent that the cortex col¬ 
lapses under hydrostatic forces and represents a convenient 
method 132^341 to account for osmotic regulation present 
in cells ED- Cortex elements are assumed to follow over¬ 
damped dynamics lic^a = where fic (kgm“^s“^) 

is the cortical drag coefficient. 

We first consider an impermeable elastic membrane that 
moves with the fiuid velocity satisfying a no-slip boundary 
condition. The fiuid velocity u and the pressure p satisfy 
Stokes equation with the incompressibility constraint 

/iAu = Vp — f, V • u = 0, (2) 

where p is the fiuid viscosity and f is the body force per 
unit volume (Nm“^) that can be calculated by spreading the 
force density from solid (Lagrangian r^^) to fiuid (Eu¬ 
ler r) coordinates as f(r) = where 

Fmfe = IS the two dimensional dis¬ 

cretized delta function. Eq. Q is solved using the regular¬ 
ized Stokeslet method ED as in Ref. E3 After the velocity 
has been calculated in Euler coordinates, we evolve the mem¬ 
brane nodes as We take special care to correct 

to enforce volume conservation, a common problem 
of the immersed boundary method I^I^IJTII . 


Cortex contraction and healing are the main driving forces 
of the blebbing activity. In our model, we assume that the cor¬ 
tex is pre-stretched which we impose by choosing a value for 
the equilibrium distance between the nodes Cc that is smaller 
than their initial distances. Bleb nucleation in our model oc¬ 
curs stochastically due to the randomness of the stiffness of 
the bonds representing the membrane-cortex interface. Lo¬ 
cal detachment is then implemented by setting to zero the 
bond stiffness kmc if its stretching {rma \ is above a threshold 
that we set at 0.1/. Einally, to implement interface healing, 
we assume that each cortex node i, associated to a removed 
interface bond, moves towards the membrane with velocity 
V = t'crmci/kmcj Until it attaches again when |ri„cj < I 
IZ7l . When disconnected nodes become connected again, we 
assign new random values to the stiffness kmc of the bond. We 
observe that the strength of disorder (i.e. the wideness of its 
distribution) controls bleb nucleation which in turn may cause 
the presence of a large number of blebs, occurring simulta¬ 
neously. As in previous models 125^281 . we neglect addi¬ 
tional time-dependent effects due to viscoelasticity and actin 
turnover in the cortex 1361 . 

We introduce permeability into the model by following Ref. 
l29l . We assume Darcy Law and impose a porous slip velocity 
normal to the membrane given by 


, , KdP K[P] 

=-^ ~- 

p on pa 


(3) 


where K is permeability, p is the viscosity, a is membrane 
thickness and [P] is the pressure jump. Using the normal 
stress jump condition [P] = the porous velocity 

Up reads i4p(rmj = -aF • n/AS'(rmJ where a = K/{pa) 
(m^skg“^) and F is the force on the node. Einally, we reach 
the equation used in the simulations Vmk = u^^^^{rmk) + 

Up(rmJ. 

Simulations: To simulate the model, we assume a square 
fiuid domain that we discretize using a square grid with a dis¬ 
cretization step dx = L/Ae and dy = L/Ae, where L is 
the length of the domain in each direction and Ae x Ae is 
the number of Eulerian coordinates. The fiuid domain covers 
both the inside and outside of the membrane and has an area 
of 100 xlOOpm. The discretization steps dx and dy are of 
size 1 pm. 

The membrane permeability in zebrafish embryos has been 
measured experimentally and is reported to be in the range 
3 X 10“^^ < a < 2 X 10“^^m^s/kg depending on the de¬ 
velopmental stage 1^ . Here, we perform numerical simu¬ 
lations under different values of a, ranging from a = 0 to 
Of = 8 • 10“^^m^s/kg, to account for AQP overexpression 
and knockdown (see Table SI for a complete list of parame¬ 
ters). When a is in the physiological range, we observe re¬ 
alistic bleb formation and retraction (Eig. and Movie S5). 
When we delete membrane permeability, setting a = 0 and 
enforcing strict cell volume conservation, bleb activity is sup¬ 
pressed (Eig. and Movie S6). We can relate the simula¬ 
tions to the experimental results by noticing that the cell area 
A, the two dimensional analogue of the three dimensional cell 
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FIG. 3: (Color online) Numerical simulations show that membrane 
porosity is needed for blebbing. a) A series of snapshots of the nu¬ 
merical simulations of bleb formation and retraction. The color rep¬ 
resents the local fluid pressure: red for positive and blue for negative 
pressures while the cell membrane is green, a) Results obtained with 
physiological permeability: a = 2 • 10“^^m^s/kg. b) For a = 0, 
bleb formation is impaired, c) The cumulative distribution of rela¬ 
tive area fluctuations for different permeabilities, d) The standard 
deviation of the distribution of relative cell areas as a function of 
permeability. The physiological range is depicted in grey. 


volume V, fluctuates more or less when the permeability is in¬ 
creased or reduced (Fig. [^d), in correspondence with AQPs 
over-expression or knock-down, respectively ll23ll . Both ex¬ 
periments and simulations suggest that blebs occur as long as 
the fluid is able to flow sufficiently fast through the membrane. 

Discussion: Our model allows us to better understand why 
volume fluctuations are crucial for blebbing. Acto-myosin 
driven cortex contraction leads to a shrinkage of the cell by 
squeezing some water outside. Indeed the fluid pressure in¬ 
side the cell is initially larger than the one outside (see Fig. 
[^). The contraction of the cortex induces stretching in the 
membrane-cortex linkers leading the membrane to buckle 
1391 . Buckling provides an effective way for membranes to 
avoid considerable elastic compression and is associated with 
a structural softening of the system |40ll4Tl. Furthermore, it 
allows to generate large membrane deflections needed to form 
a bleb. In differentiated cells additional membrane surface 
can be obtained by disassembling caveolae l42l . but this can 
not happen in PGCs where caveolin is not expressed. When 
the interface fractures, the mechanical stress on the detached 
part of the membrane decreases, but it increases on the inter¬ 
face that is still attached, inducing crack propagation. As the 
bleb expands, the fluid pressure inside the cell is reduced (see 
Fig. [^) leading to an inflow of water. In Fig. |^, we display 
the spatio-temporal evolution of the pressure jump across the 
membrane, showing large fluctuations in correspondence to 
bleb formation 1221 . These pressure spikes, whose distribu¬ 
tion is long-tailed (Fig. |^), are a manifestation of the stress 
concentrations around cracks and are needed to account for 



FIG. 4: (Color online) a) Spatio-temporal evolution of the pressure 
jump across the membrane computed from the normal stress jump 
condition in numerical simulations with a — 2- 10“^^m^s/kg. The 
pressure drop is strongly enhanced in localized regions correspond¬ 
ing to blebs, b) The corresponding distribution of pressure jumps 
displays long tails. 


the observed volume fluctuations, because otherwise the av¬ 
erage pressure jump generated by a uniform cortex contrac¬ 
tion (around 10^ Pa) would not displace sufficient amount of 
fluid during the short lifetime of a bleb (141 [35l. Healing of 
the membrane-cortex interface eventually leads to bleb retrac¬ 
tion and to an increased fluid pressure inside the cell. The 
mechanism described above does not work for an imperme¬ 
able membrane: isochoric buckling is possible in principle 
but, in addition to bending, it necessarily causes considerable 
stretching which is energetically expensive (3^ . Thus, when 
the interface fractures the bleb does not form and interface 
delamination takes place without localized membrane expan¬ 
sion. 

Our numerical results show that volume fluctuations during 
blebbing are due to mechanically induced highly non-uniform 
pressure drops across the membrane. The same mechanism 
could also explain the experimental observation that blebs are 
nucleated preferentially in regions of negative membrane cur¬ 
vature da. Changes in osmotic pressure gradients could also 
contribute to the process, as suggested previously MMM, 
but are not explicitly included in our model. Future experi¬ 
mental work will clarify if this and other assumptions present 
in our model are correct, but our results should stimulate both 
new experiments and the development of more elaborate the¬ 
ories and models. This would also allow to better understand 
the role of transmembrane water transport for other cellu¬ 
lar protrusions, given that past experimental results relate the 
presence of aquaporins to the formation of lamellipodia El 
and filopodia l20l . The present methodology provides the ha- 
sis for a physical explanation of this broad class of phenom¬ 
ena. 
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SUPPLEMENTARY INFORMATIONS 
Zebrafish maintenance 

Zebrafish (Danio Rerio) of the AB and AB/TL genetic 
background were maintained, raised and staged as previously 
described 14511461 . 


Constructs and mRNA synthesis 

For protein over-expression in germ cells, the mRNA was 
injected into the yolk at one-cell stage. Capped sense RNA 
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was synthesized with the mMessage mMachine kit (Ambion, 
http://www.ambion.coni/index.html). To direct protein ex¬ 
pression to PGCs, the corresponding open reading frames 
(ORFs) were fused upstream to the 3UTR of the nanosl 
(nosl-3UTR) gene, facilitating translation and stabilization of 
the RNA in these cells i?7l . For global protein expression, 
the respective ORFs were cloned into the pSP64TS ector that 
contains the 5 and 3 UTRs of the Xenopus Globin gene. The 
injected RNA amounts are as provided below. The following 
constructs were used: 

- YPet-YPet-RasCAAX-nos-1 (240 pg.) was used to la¬ 
bel membrane in germ cells. 

- Lifeact-pRuby-nos-1 (240 pg.) was used to label actin 
in germ cells. 

- DN-ROK-nos-1 (300 pg.) was used to interfere with 
ROK function in PGCs 

- Aqpla-nos-1 (300 pg.) was used to over-express 
aquaporin-la in PGCs 

- Aqp3a-nos-l (300 pg.) was used to over-express 
aquaporin-3a in PGCs 

- AqplaEGFP-nos-1 (360 pg.) was used to visualize the 
subcellular localization of aquaporinla in PGCs 

- Aqp3aEGFP-nos-l (300pg.) was used to visualize the 
subcellular localization of aquaporin3a in PGCs 

Morpholino knockdown 

The morpholinos for knocking down protein translation 
were obtained from GeneTools, EEC http://www.gene- 
tools.com/ The following sequences were used: Aquapoinla: 
5 AAGCCTTGCTCTTCAGCTCGTTCAT3 (injected at 
400/iM); Aquaporin 3a: 5 ACGCTTTTCTGCCAACC- 

CATCTTTC 3 (injected at 400/iM);. Eor the control, standard 
morpholino 5CCTCTTACCTCAGTTACAATTTATA 3 was 
used. 


Live Imaging of germ cells in zebrafish embryos 

Time-lapse movies of blebbing cells in live zebrafish em¬ 
bryos were acquired with the Zeiss LSM710 bi-photon micro¬ 
scope using one-photon mode. The 20x water-dipping objec¬ 
tive with the numerical aperture 1.0 was used. The bit depth 
used was 16 and the scanning speed ranged between 150 to 
250 ms/frame for fast imaging of bleb formation. 

Image processing 

Images were preprocessed with Eiji software to elimi¬ 
nate the background. The Bleach Correction tool (EMBL) 


was used to correct for the reduction in fiuorescence in¬ 
tensity during prolonged time-lapse movies. Sequences of 
of image stacks were then processed using the 3D Ac¬ 
tive Meshes algorithm |[24ll implemented in the Icy soft¬ 
ware http://icy.bioimageanalysis.org/ The algorithm per¬ 
forms three-dimensional segmentation and tracking using a 
triangular mesh that is optimized using the original signal 
as a target. Erom the resulting three dimensional mesh one 
can then measure the cell volume its surface area. Three di¬ 
mensional rendering of the meshes was done using Paraview 
(http ://www.paraview.org/). 

Statistical analysis 

Statistical significance was evaluated using Kolmogorov- 
Smirnov tests implemented in custom made python codes. 

Test of the accuracy of cell volume evaluations 

To confirm that the volume fiuctuations observed in the ze¬ 
brafish cells in vivo are real, we have to exclude possible sys¬ 
tematic errors induced by the three dimensional mesh recon¬ 
struction algorithm 1241 . Our experimental analysis shows 
smaller volume fiuctuations for cells where the bleb forma¬ 
tion has been suppressed (DN-ROK and AQP- mutants), with 
respect to wild type cells (WT) or those for which blebbing 
has been enhanced (AQP-f). In particular, as shown in Pig. 
1(a), WT and AQP-f cells display a more complex morphol¬ 
ogy when compared to DN-ROK or AQP-, which appear in¬ 
stead to have a rounded shape. Hence, the first question to 
be answered is whether in presence of complex morphologi¬ 
cal shapes, the algorithm introduces a systematic bias in the 
measured cell volume. In other words the question is: does 
the algorithm produce a larger error while calculating the vol¬ 
ume of blebbing cells, with respect to those where blebs are 
absent? 

To answer to this question we generate a set of synthetic 
ellipsoidal cells whose volume is in the range of the zebrafish 
cells analyzed in our experiments (^ 2000 — 3000/im^), 
as shown in Pigj^ Synthetic cells are generated through 
the Image! software (http://imagej.nih.gov/ij/) by creating 3D 
stacks having the shape of ellipsoids with semi-axes ay 
and az. Initially, we set the voxel size to 0.13/im along the 
three directions (see PigjS^a)). The volume is then calcu¬ 
lated both according to the formula Veil = /‘^a^ayaz and 

by counting the number of voxels belonging to the synthetic 
cell, Vwp^ The relative error between the two estimates is 
^ _Q SO that we can safely consider the 
estimate as the real volume of the synthetic stacks gener¬ 
ated. The main sources of error in analyzing confocal image 
stacks stems from the anisotropic voxel. In our experiments 
the resolution along the z direction is 0.78/im, while it is 
0.13/im in the xy plane. To reproduce the voxel anisotropy 
in the synthetic stacks, we select just one single xy plane 
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every 6 composing the original z-stack. A resulting typical 
cell is shown in Fig|S3[b). We then extract the mesh of this 
newly obtained stack (see FigjS^c)) and calculate the en¬ 
suing volume Vno-bieb- Our set of synthetic cells consists 
of 35 ellipsoids of different semi-axes, for which we calcu¬ 
late the true volume V^p reported in Fig|^ Then, each el¬ 
lipsoid is first processed by the anysotropic voxelization in 
the z direction, and subsequently analyzed by the 3D Active 
Mesh algorithm. The volumes of the extracted meshes are 
reported in Fig|^ (Vno-bieb)- H is apparent that the algo¬ 
rithm systematically underestimates the volume by roughly 
j _ / Vno-bieb-v,^p K ^ —14%. We checked that is error is 

' ^wp ' 

greatly reduced for isotropic voxels. 

A constant systematic error is not worrying, since we are 
only interested in changes in volume and all the images have 
the same voxel anisotropy and therefore the same error. Be¬ 
fore addressing the fiuctuations of this error, we focus on pos¬ 
sible spurious changes in the measured volume induced by a 
change in shape. For each of the synthetic ellipsoidal cells, 
we create a synthetic cell with the same volume but present¬ 
ing 1, 2, or 3 blebs on the surface. Blebs are generated as 
spherical caps of different radii with centers placed randomly 
on the ellipsoid surface (see FigjS^d)). We then perform the 
same anysotropic voxelization of the original z-stack done for 
the plane ellipsoids (see FigjS^e)). From this image we ex¬ 
tract the active mesh (FigjS^f)) and calculate its volume. The 
volumes of the meshes of synthetic cells with blebs, Vueb^ are 
displayed in Figj^ One can only see a very small difference 
between the values of Vueb and Vno-bieb, but they both ap¬ 
pear underestimate the true value V^p by about 6 —14%. 

What is surprising, however, is that cells with blebs appear to 
approximate the real volume slightly better than cells without 
blebs. To confirm this, we report in Fig|S5|the relative volume 
fiuctuations as a function of the measured mesh 

y no-bleb 

surface fiuctuations ^bieb-^no-bieb errors were made by 

^no—bleb 

the algorithm in estimating the synthetic volumes one would 
expect ^bie^-Vno-bieb _ Q pail* of cells were constructed 

with the same volume but different shapes. If complex shape 
with blebs would lead to an overestimation of the volume with 
respect an ellipsoidal cell with no blebs, one would expect 
Vbieb-Vno-bieb bzeb posltivcly corTclatcd. 

To the contrary, the linear regression of the data in Fig|S5] 
shows a small but clear anti-correlation between volume and 
surface fiuctuations. This is in contrast with experimental re¬ 
sults showing that changes in shape are positively correlated 
with changes in volume (Fig. [S7| ). Hence, this result can not 
be considered an artefact of the measurement but a real feature 
of the cells. 

The previous analysis clearly demonstrates that the volume 
fiuctuations observed in our experiments do not depend on 
the shape of the cells, but the algorithm systematically un¬ 
derestimates the volume of both of about 5 —14%. The 

next question is whether these systematic fiuctuations are of 
the same order of magnitude as the observed ones. Are the 
fiuctuations reported in Fig. 1(b) real or just an artefact in¬ 


troduced by the mesh reconstruction algorithm? This ques¬ 
tion is particularly compelling in the case of WT and AQP-i- 
cells, since for DN-ROK and AQP- cells we can accept that 
the volume might remain constant. To answer to this question 
generate a set of 120 synthetic cells, 60 with blebs randomly 
placed and of different sizes, and 60 without blebs, each cell 
having its own real volume V^p calculated with Image!. We 
then process each synthetic cell according to the protocol pre¬ 
viously outlined: anysotropic voxelization in the z direction 
and subsequent mesh analysis. Finally we calculate the fiuc¬ 
tuations ^ and its cumulative distribution 

V _ y wp 


In Fig 


S6 


sponding cumulative 


we compare P with the corre- 

distributions of WT, AQP-I-, DN-ROK 
and AQP- cells, once it has been shifted by the average sys¬ 
tematic volume bias 6 —14%. This figure shows that sys¬ 

tematic errors made by Icy in the estimation of the cell vol¬ 
ume, are compatible with the observed volume fiuctuations of 
DN-ROK and AQP- cells, but not with those AQP-i- and WT 
cells, which instead appear to be significantly larger. We thus 
conclude that the difference in volume fiuctuation between 
AQP-fAVT cells and DN-ROK/AQP- is not an artefact of the 
analysis. 


Volume conservation in numerical simulation 

By performing numerical simulations, we notice that a 
straightforward formulation of the method suffers from poor 
volume conservation. This general drawback of the immersed 
boundary method has been already pointed out by many other 
authors in different contexts 12911371. To overcome this prob¬ 
lem we implement the method proposed in Ref. Ell and 
enforce the incompressibility constraint on the discrete La- 
grangian grid in the weak sense: 

N 

y]u(rmJ-n(rmJAS'(rmJ = 0, (4) 

k=l 

where n(rmfe) is the outward unit normal to the membrane at 
the position and A5'(rmfc) is a discrete measure of the ar- 
clength in the actual configuration at the position . The 
above constraint is satisfied by adding a corrective term to 
the equation of motion where = 

^(i*mfc) — satisfies the incompressibility constraint 

and M is given by 

1 ^ 

M = ^ --^u(rmj •n(rmJAS'(rmJ. (5) 

Ei=iA5(r„,J 

With this correction, we observe that the incompressibility 
constraint is satisfied and the cell volume is perfectly con¬ 
served. 


SUPPLEMENTAL FIGURES 
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Figure SI: AQP3 is expressed by PGCs. A PGC expressing EGFP 
fusion of Aquaporin-3a and an RFP-tagged membrane marker. Scale 
bar is 10/im. 



Figure S2: Aquaporin knockdown and overexpression does not in¬ 
duce significant changes in the average cellular volume. A similar 
result holds for the DN-ROK mutant. The results show the disper¬ 
sion of the data, the average and the standard error. Statistical signif¬ 
icance (p values) is evaluated according to the Kolmogorov-Smirnov 
test. 


Figure S3: Two example of synthetic image stacks representing cells 
without (a-c) and with (d-f) blebs. The original stacks of equal vol¬ 
ume (a and d), are first transformed removing a set of planes (1 every 
6) to obtain an anisotropic voxel corresponding to the experimental 
resolution (b and e), i.e. 0.13/im in the x,y directions and 0.78/im 
along z. The resulting stacks (b and e) are analyzed to obtain a three- 
dimensional mesh (shown in c and f) 



Figure S4: Volume comparison of synthetic cells. We compare 
the measured volumes of cells without blebs (Vno-bieb) with those 
corresponding to cells with blebs (Vbieb) and with the true volumes 
(Ywp), which is the same for each pair of cells. The results show a 
large constant systematic error 6 ~ —14%, and small fluctuations 
for cells with and without blebs. 
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Figure S5: The relative difference in the measured volume for pairs 
of cells of the same true volume but different surface due to the pres¬ 
ence or absence of blebs. The data show fluctuations of 1% and a 
small negative correlation between volume and surface changes. 


Figure S7: Volume and surface fluctuations are correlated. Principal 
component analysis of volume and surface relative values. Scatter 
plots for a) WT, b) DN-ROK, c) AQP-i- and d) AQP- are reported 
together with an ellipse with axis given by the two eigenvectors of 
the cross-correlation matrix, whose amplitude is reported in panel e) 
and f) for the largest and smaller axis. The dashed line represents 
the expected result for the ideal case of the isotropic deformation of 
a sphere. 



Figure S6: The cumulative distribution of relative error fluctuations 
for the volume of a large number of synthetic cells with and without 
blebs compared with experimental measurements. The cumulative 
distribution of synthetic cells is shifted by —(5 to allow a visual com¬ 
parison with the experimental quantities. V corresponds the time av¬ 
eraged volume for AQP-I-, AQP-, DN-ROK and WT cells, whilst it is 
V = Vwp for synthetic cells. Volume fluctuations for WT and AQP-i- 
cells are signiflcantly larger than those observed in synthetic cells, 
whereas DN-ROK and AQP- mutant cells volume fluctuations seem 
to be compatible with the systematic errors induced by the mesh al¬ 
gorithm. 



Figure S8: Time evolution of the pressure drop across the mem¬ 
brane from numerical simulations.A representative example of the 
evolution of the pressure drop in numerical simulations shows that 
the average is very different from the maximum (top). Furthermore, 
the standard deviation of the distribution fluctuates intermittently in 
time in correspondence to the blebbing activity (bottom). Results are 
obtained for a permeability a = 2 - 10“^^m^s/kg. 
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Figure S9: A schematic representation of the bleb formation pro¬ 
cess. a) The cortex contracts squeezing water outside of the cell, b) 
The membrane buckles and the cortex-membrane interface fractures, 
c) The bleb expands as the interface fails and water flows inside the 
cell as the internal fluid pressure is relieved. 
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SUPPLEMENTAL TABLES 



WT DNROK AQP- AQP+ 

p-value (AV/V) 
p-value (AS/S) 

0.11 0.21 0.17 0.99 
0.68 0.75 0.7 0.13 


Table S1: Results of statistical significance tests for validity of Gaus¬ 
sian statistics for volume and surface fiuctuations. We report the p- 
values obtained from the Kolmogorov-Smirnov test. A small p-value 
(e.g. p < 0.01) would imply that we can reject the hypothesis that the 
distribution is described by Gaussian statistics. In the present case, 
the p-value is large indicating that a Guassian distribution provides a 
good fit to the data. 



WT / WT/ WT/ AQP-/ AQP-h/ AQF+/ 

DNROK AQP- AQF+ DNROK AQP- DNROK 

p-value (AV/V) 
p-value (AE/E) 

5 • 10“'’ 4 • 10“^ 0.25 062 8 • 10“'’ 8 • 10“'’ 

0.001 0.008 0.09 0.98 6 • 10“® 6 • 10“® 


Table S2: Results of statistical significance tests for the compari¬ 
son between WT, DNROK,AQP-i-,AQP- cells in the case of volume 
and surface distributions. We report the p-values obtained from the 
Kolmogorov-Smirnov test. A low p-value (e.g. p < 0.01) indicates 
that we can reject the hypothesis that the two data sets are described 
by the same distribution. 


Symbol 

Quantity 

Value 

Reference 

Tm 

membrane radius 

25 /im 

inr 

Tc 

cortex radius 

24 pm 


€c — €m 

rest length of cortex/membrane elements 

1.3 pm 

— 

A 

regularization parameter 

0.5 X e 

EH 

km 

membrane stiffness coefficient 

6 X 10“® Nm“^ 

O 

kc 

cortex stiffness coefficient 

9 X 10“® 

O 

Bm 

membrane flexural rigidity 

4 X 10“^° J 

ca 

Be 

cortex fiexural rigidity 

2.8 X 10“^® J 

ca 

kmc 

cortex-membrane interface stiffness 

25 X 10“® Nm“® 

ca 

pic 

cortex drag coefficient 

10“^ kgm“^s“^ 

Eli 

l^c 

cortex healing speed 

6 X lO-^ ms“^ 

l27l 


Cytosolic viscosity 

10“^ kgm“^s“^ 

l27l 

a 

membrane permeability 10“^ 

- 10“^ X 10“^2 m®skg- 

ED 


Table S3: Parameters employed in numerical simulations 
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SUPPLEMENTAL MOVIE CAPTIONS 


PGC in AQP- conditions. The movie is obtained using 
an average intensity 3D projection in imageJ. 


Movie SI A representative example of the time evolution ^ 

j v j • Movie S4 A simulation ol the computational model using a porous 

PGC m WT conditions. The movie is obtained using ^ ^ ^ ^ ^ 

• ^ • T membrane. The color represents fluid pressure (see Fig. 

an average intensity 3D projection m imageJ. ^ r v & 

3a}. 

Movie S2 A representative example of the time evolution of a 

PGC in DNROK conditions. The movie is obtained us -^ovie S5 A simulation of the computational model using an im- 
ing an average intensity 3D projection in imageJ. permeable. The color represents fluid pressure (see Fig. 

3b). 

Movie S3 A representative example of the time evolution of a 



